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Abstract 

Spider dragline silk is one of the toughest materials known and understanding the hierarchical 
structure is a critical component in the efforts to connect structure to function. In this paper, 
we take the first step in elucidating the hierarchical fractal structure of /I-sheet nano-crystallites, 
which form a robust self-similar network exhibiting an non-linear mechanical property. A combined 
small angle X-ray scattering (SAXS) and wide-angle X-ray scattering (WAXS) study of the nano¬ 
crystalline component in dragline silk fibers from several species of spiders including, Latrodectus 
hesperus, Nephila clavipes, Argiope aurantia and Araneus gemmoides is presented. SAXS structure 
factors exhibit a ‘lamellar peak’ in the q-range from 0.60 to 0.82 nm"^ for various spider dragline 
silk fibers, indicating the presence of strong nano-crystal ordering on the >10 nm length-scale. The 
stochastically reconstructed electron density maps indicate that the /3-sheet crystals are hierarchi¬ 
cally structured as mass fractals and that nano-crystals tend to form 10 to 50 nm sized clusters 
with long-range crystalline ordering. This nano-crystal ordering along the fiber axis also helps 
to explain the difference between axial and radial sound velocities recently measured by Brillouin 
spectroscopy. 
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I. INTRODUCTION 


Spider dragline silk fibers have excellent reversible extensibility and high tensile strength 
Dia. Extensive studies on the molecular structure of spider silk have been performed over 
the years and it’s widely believed that oriented anti-parallel /5-sheet nano-crystals are the 
key contributor to spider silk’s excellent mechanical properties [3] . The molecular structure 
and chemical composition of spider silks have been studied extensively by solid state nuclear 
magnetic resonance spectroscopy and these studies have shown that a large fraction of the 
amino acid sequences are poly-(Gly-Ala) and poly-Ala repeats [IHTj, which form rigid anti¬ 
parallel /5-sheet nano-crystals through the periodic hydrogen bond assemblies [3[|HIE]. The 
less ordered amino acids are in the form of random-coil like helical secondary structures in 
which the /5-sheet nano-crystals are embedded. The elastic and random-coil like a- and 3io- 
helical structures are abundant and occupy a large volume of the fiber body, acting as the 
interconnections among the rigid crystals HDl. The crystal structure and physical size of 
individual /5-sheet nano-crystal has also been resolved by X-ray diffraction studies. Past 
WAXS studies, as well as this work, have confirmed that typical /5-sheet crystals have an 
orthorhombic unit cell with their physical sizes ranging from 2 to 4 nm, when produced at 
the natural extrusion speed mm- 

While the nano-scale dimensions of the /5-sheet crystals have been well studied, their 
hierarchical structures and relation to the macroscopic mechanical properties still hasn’t 
been solved. Being the most rigid objects in the spider dragline silk, the /5-sheet crystals 
play a crucial role in determining its mechanical properties. If we assume the /5-sheet 
crystals to be the building blocks of the cylindrical-shape fiber, then its physical size, inter¬ 
crystal distances and long-range packing pattern will be the key parameters that define 
the macroscopic mechanical properties. Fortunately, several microscopy studies have gained 
insight on the crystallite structures of the spider silk fibers. Scanning electron microscopy 
(SEM) studies have shown that the texture of ion-etched silk fiber’s is rather rough, as 
scattered crystalline-rich regions about the size of 20 nm to 50 nm have been observed 
across the silk fibril ra- Transmission electron microscopy (TEM) image shows relatively 
large crystallites on the scale of 70 to 120 nm are embedded in the amorphous matrix m- 
To date, there is still no consistent model to describe the hierarchical structure of these 
large crystalline regions, as the dominant scattering centers, i.e. /5-sheet crystals, only span 
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several nanometers in all three dimensions. By analyzing of the SAXS structure factor 
using a stochastic reconstruction method, we provide evidence that the crystalline structure 
of spider dragline silk hber is mass fractal, accompanied by dense clustered packing of the 
nano-crystals. The ‘large crystals’ that span up to 70 nm are composed of highly oriented 
and closely interlinked /3-sheet crystals. In this study, we will present a combined WAXS and 
SAXS analysis of the crystalline phase in spider silks followed by the modeling of the spacial 
packing of the /3-sheet crystals that reveals the underlying morphology of the crystalline 
structure within the spider dragline silk fibers. 


TABLE I. Lattice parameters, nano-crystal sizes and inter-crystallite distance. 

Lattice parameters From (200) From (120) lamellar peak d-spacing 


Fiber 

“(A) 

6(A) 

c(A) 

a 

29 

Ti(A) 

a 

29 

7-2 (A) 

d(A) 

din 

dlT2 

L. hesperus 

10.46 

9.62 

6.88 

1.102 

15.18 

21.3 

0.824 

18.24 

26.4 

77.1 

3.6 

2.9 

N. elavipes 

10.70 

9.73 

6.86 

1.073 

15.21 

20.3 

0.960 

18.26 

22.7 

83.2 

4.1 

3.6 

A. aurantia 

10.52 

9.74 

6.92 

1.107 

15.06 

19.6 

0.933 

18.13 

23.4 

108.7 

5.5 

4.6 

A. gemmoides 

10.58 

9.67 

7.09 

1.023 

15.09 

21.3 

0.902 

18.10 

23.7 

106.6 

5.0 

4.5 


II. RESULTS AND DISCUSSION 


Fig.l shows the WAXS and azimuthal integrated SAXS prohle of L. hesperus (Black 
Widow) dragline silk hber. The diffraction pattern is divided into two distinct regions: the 
center small-angle region (q<1.3 nm'^) and the outer wide-angle region. As shown in Fig.la, 
the diffraction spots have been assigned with Miller indices from which we have calculated 
the unit cell parameters basing on an orthorhombic unit cell model [12]. To retrieve the 
crystal sizes from X-ray data, we integrated the wedge shape equatorial region containing 
rehections (200) and (120) (Figj^), and then applied a 1-D Gaussian peak htting procedure 
[see Fig.SI]. The full maximum at half width (FMHW) of the htted Gaussian peaks are 
evaluated using the Scherrer equation na 


K\ 


T = 


( 1 ) 


(3 cos 6 

where K = 0.9, A = 1.38 A, (3 is the FWHM of htted peak and 9 is the Bragg angle, to 
calculate the nano-crystal physical dimension r. The orthorhombic unit cell and crystal 
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FIG. 1. A combined small angle X-ray scattering (SAXS) and wide-angle X-ray scat¬ 
tering (WAXS) image, (a) WAXS pattern of L. hesperus (Black Widow) major ampullate 
(dragline) silk. The wide-angle diffraction spots have been assigned with Miller indices, from 
which the unit cell parameters have been calculated using an orthorhombic crystal model. Gaus¬ 
sian peak fitting was applied to the wedge shaped region (white line) and the Scherrer equation 
was used to determine the average nano-crystallite dimensions. The center region (q<1.3 nm'^) 
corresponds to the SAXS scattering pattern, (b-d) Samples of N. clavipes, A. aurantia and A. 
gemmoides show similar wide-angle scattering patterns, (e) Azimuthal integration of scattering 
intensity from L. hesperus dragline silk fibers. The lamellar peak is located at q=0.82 ±0.01 nm'^. 
The inset shows the SAXS structure factor S{q) on log-log scale, where the ‘matrix knee’ represents 
the intermediate length scale (1 nm to 200 nm) and exhibits linearity with a slope of 2.16. 
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sizes are summarized in Tablej^ The crystal size in all three dimensions (a, b, c) shows very 
small variations for the four species of dragline silk examined here. The crystal sizes have 
a narrow distribution, namely from 19.6 to 21.3 A for that calculated from (200) reflection 
and from 22.7 to 26.4 A for the (120) reflection. Although different species of dragline silk 
shows signihcant differences in alanine content and mechanical propertiesjl], their building 
blocks are surprisingly similar in terms of physical appearance. Fig.le shows the small- 
angle scattering prohle of L. hesperus (Black Widow) major ampullate (dragline) silk hbers, 
where the characteristic lamellar peak is located at 0.82 nm'^. The characteristic correlation 
peaks between 0.5 and 1.2 nm'^ have been observed frequently in polymer and copolymer 
materials dS]. The correlation peaks manifest certain long-range ordering of the crystalline 
phases and such ordered state is crucial to the functional behavior of the material, such as 
water channels morphology observed in nahon polymers HD. Therefore, we propose that 
it’s the intermediate length-scale morphology of the /5-sheet crystals that determines the 
macroscopic properties such as mechanical strength, elasticity and thermal conductivity. 



FIG. 2. Correlation of lamellar peak positions to Alanine content. Alanine in spider silks 
is shown to primarily occur in /3-sheet nano-crystals and an increase in alanine content in the 
spider silk protein is a contributing factor to increased crystallinity associated with the nanoscale 
clusters which give rise to the lamellar peaks. The alanine content data was previously retrieved 
from solid-state nuclear magnetic resonance experiments. 

Since /3-sheet crystals are identihed as the only kind of strong scattering centers in spider 
dragline silks, it’s natural to assume /3-sheet crystals as the origin of the SAXS lamellar 
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peaks. The packing pattern of the crystals has unique intermediate-range ordering, which 
leads to the appearance of strong lamellar peaks in the g < 1 nm'^ range of the SAXS 
structure factor [TSl [12]. For strong correlation peaks to appear, it’s important that the 
scattering centers maintain a closest approach distance Rca between each pair, as has been 
discussed by Yarusso and Copper dg. In real space, the closest approach distance Rca 
limits the average spacing between each adjacent pair of scattering centers that reside in 
the amorphous backbone of biopolymer. Reflected in reciprocal space, the combination 
of closest approach distance and long range ordering will generate correlation peak in the 
small-angle scattering regime. For different species of spider dragline silks examined here, 
all of them exhibit correlation peaks in the range of 0.6 to 0.9 nm'^. This indicates that 
the /3-sheet crystals are not adjacent to each other statistically, but rather maintain an 
average closest distance among them, which can be quantihed by the characteristic of SAXS 
lamellar peak. The SAXS signal exhibits an elliptical streak, elongated along the meridian 
direction. This indicates that the /3-sheet nano-crystals have a rectangular shape with long 
axis parallel with hber axis, as illustrated by Fig.S2 in supplementary material. Combining 
these information, we constructed a 2-dimensional electron density map in which the beta- 
sheet crystals are initialized with constrained orthorhombic geometry. The d-spacings of the 
lamellar peaks were taken as the average inter-crystal distance in constructing the initial 
model. This in turns determines the density of the /3-sheet crystals in each silk species. The 
linear correlation between the lamellar peak position and alanine content further strengthens 
this assumption. 

Past solid state NMR experiments have shown that the fraction of alanine content, which 
primarily occur in the /3-sheets, vary signihcantly among different species HEDI. By cross 
comparing the data from x-ray scattering and NMR, we have found the alanine content 
to be linearly correlated to lamellar peak’s position in g-space, as shown in Figj^ With a 
relatively small inter-crystal distance, such as 83.2 and 77.1 A measured from N. clavipes 
and L. hesperus respectively, we hnd that the /3-sheet crystals are more densely packed and 
thus leading to a higher fraction of beta-sheet crystals than the other two samples, namely 
A. aurantia and A. gemmoides (Fig.S3). 

The experimental and numerically simulated SAXS structure factors S{q) are shown in 
Figj^. The lamellar peaks are present in all hber samples and they range from 0.6 to 
0.9 nm'^ with variations in the peak intensity. The lamellar peaks of A. aurantia and A. 
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FIG. 3. Structure factors S{q) and pair correlation function from reconstructed struc¬ 
ture. (a) The numerically simulated structure factor S{q) (—) fits the experimental SAXS S{q) 
(---). The small-angle lamellar peaks’ positions, from L. hesperus (top) to A. gemmoides (bot¬ 
tom), are correspondingly located at 0.83, 0.78, 0.60 and 0.61 nm'^ with errors of approximately 
TO.Olnm'^. The low-g' region of the structure factor S{q), i.e. the matrix knees, exhibit linearity 
in all cases and the slopes range from -2.1 to -2.7 on the log-log scale, (b) Pair correlation function 
P(r) calculated from the reconstructed electron density maps. The correlation function P(r) were 
calculated from a population of around 5000 crystals and the curves were numerically smoothed. 
The intermediate range crystalline ordering is reflected as the multiple correlation peaks observed 
on the P(r) curves in the range of 7 to 40 nm. 

gemmoides are relatively close in q-space with a minuscule difference Ag=0.01 nm'^, while 
L. hesperus and N. clavipes are slightly further apart with a difference of Ag=0.05 nm'^. 
The structure factor curves are very well reconstructed across the entire collected q-range 
by the stimulated annealing reconstruction method. The lower bound of the S{q) is limited 
by the detector coverage on 14-ID-B beam line while the upper bound is cut at q=1.2 nm'^, 
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FIG. 4. Time evolution of the structure factor and the pair-correlation function P(r) 
during the stochastic reconstruction on silk species A. gemmoides . (a) From bottom 
to top, as the stimulated annealing temperature T drops, the calculated structure factor (—) 
converges to the experimental structure factor (---), reducing the pseudo-energy E as defined in 
Eq.4. (b) Initially aX E = 0.4256, the P{r) function is absent from any correlation peak. As 
the stimulated annealing algorithm proceeded, the model build up intermediate range crystalline 
ordering which was reflected by the correlation peaks at 12, 15, 21, 26, 32, 37 and 42 nm on the 
top curve {E = 0.0021) . The correlation function P(r) is sampled from a model contains 4721 
/3-sheet crystals. 


beyond which point the strnctnre factor S{q) begins to exhibit WAXS featnre. Figj^ shows 
the pair correlation fnnctions P(r) calcnlated from the reconstrncted electron density maps. 
The inter-molecular /3-sheet pair correlation function P(r) [211 E2] is dehned as 


P{r) 


2nrpoN {wf 




( 2 ) 


,where 27rrpo is the density of states in 2-D isotropic, homogeneous system, Wi is the weight¬ 
ing factor for the corresponding scattering center, and N is the /3-sheet crystal population. 
The intermediate range crystalline ordering is reflected as the multiple correlation peaks ob- 

















FIG. 5. Reconstructed electron density map of silk species A. gemmoides. (a-c) Coarse 
grained electron density map from high to low degree of coarse graining. Each map shows the 
magnified upper right quater of the previous. The texture on these three maps with drastically 
different length scales show self-similarity property. The clustering effect of /S-sheets is visually 
emphasized by the boundaries in (c). (d) The schematic illustration of the hierarchical mass fractal 
structure of the /3-sheet clusters in silk fiber, from SEM image to the basic /3-sheets network. The 
micro-fibril (100-200 nm) is composed of functional mechanical units, and each unit (40-60 nm) is 
composed of multiple crystallites-rich MaSpl proteins (10-20 nm). Within a MaSpl, it contains 
multiple /3-sheet nano-crystals (2-3 nm). Such fractal structure is mechanical robust and exhibits 
non-linear force-extension behavior. 


served on the P(r) curves in the range of 7 to 40 nm. The hrst correlation peaks appear near 
6 nm and these peaks arise from the dominant closest pair-pair interaction of the /3-sheet 
crystals. The second peaks are in the range of 7 to 10 nm, which arise from the intra-cluster 
interaction. As shown in Figj^, the second peaks is shifting right top-down, which follows 
the density and inter-crystal spacing constrain, though the exclusion region was reduced 
during the stochastic reconstruction to allow higher mobility (Fig.4). The intensity and 
sharpness of the lamellar peak reflects strength of pair correlations directly. While the A. 
gemmoides and N. clavipes have the stronger lamellar peak (Figj^), they also exhibit more 
correlation peaks in the > 10 nm range of their P{r) function. This indicates that a larger 
fraction of the /3-sheets are in long range ordered state for these two silk species. 

Figj^ shows the evolution of the stimulated structure factor and the pair-correlation 
function calculated from the reconstructed crystalline model. The correlation function shows 
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no feature in the early stage of the reconstruction, mainly due to the randomized landscape 
of the crystals and the relatively small population, which varies between 4000 to 7000. As the 
structure factor converges to the experimental value, the correlation peaks gradually build 
up, notably for the peaks at 10, 21, 31 and 42 nm, indicated by the dashed line on Figj^. 
On one hand, this dynamic proves that one can reconstruct the pair correlation function 
through reciprocal-space reconstruction. On the other hand, it consolidate the belief that 
the SAXS lamellar peaks, i.e. the correlation peaks between the q range of 0.6 to 0.8 nm'^, 
arise from the intermediate length scale ordering of the /3-sheet crystals. 

The coarse-grained electron density maps are shown in Fig|^a-c). The lighter area indi¬ 
cates the presence of a higher density of /3-sheet crystals in that region whereas the darker 
area represents the amorphous backbone, which represents diffuse X-ray scattering. The 
/3-sheet crystal distributions are initialized with a lamellar modulation when building the 
initial model. The modulation, which typically has a strip size of Niameiiar = 128 equaling to 
a physical size of 30 - 50 nm, is designated to represent the micro-hbril structure observed 
from both SEM 1131 and AFM experiments [22] • The existence of the parallel lamellar 
nano-hbrils could explain the discrepancy between the spider silks’ axial and radial sound 
velocities reported by Koski et. al [21]. Along the axial direction of hber, these parallel 
lamellar hbrils contain a high density of /3-sheet crystals with crystalline ordering up to 30 
nm, which can be observed on Figj^. The highly ordered and continuous lamellar back¬ 
bones act like phonon highways and thus support fast propagation of both longitudinal and 
transversal acoustic waves in the spider dragline silks, as observed in the recent Brillouin 
scattering measurements [21]. On the other hand, due to the modulation of the lamellar 
structure, the crystalline ordering vanishes and the /3-sheet crystal structure has a zero- 
density discontinuity along the radial direction. Consequently, the scattering of phonons is 
stronger and the measured velocities of the acoustic wave are much lower than that along 
the axial direction. 

The other prominent structural feature is the clustered packing of the /3-sheet crys¬ 
tals. During the reconstruction process, /3-sheet crystals self assemble to form high density 
crystalline-rich islands. The clustering effect exist in all length scale examined here and the 
cluster sizes change from 30 nm in Figj^ to 50 nm in Figj^, and to 100 nm in Figj^, 
increasing exponentially in accordance with the coarse graining. The reconstructed electron 
density map shows remarkable self-similarity property and is scale invariance statistically. 
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The formation of these clusters is driven by both the characteristic of the matrix knees and 
the presence of the lamellar peaks. The SAXS structure factor has power law relation with 
respect to the scattering vector q 


S{q) oc q 


(3) 


and —r is the slope of the matrix knee in the log-log plot of S{q) (Figj^). The parameter 
r is the fractal dimension in system exhibiting self-similarity [EIESIES]. For the silk hbers 
examined in this study, the r is between 2.16 and 2.68, which means the crystalline structure 
has a mass fractal (2 < r < 3) |271I2B]- Crystal clustering is the dominant form of mass fractal 
and therefore the clustering effect can be observed at drastically different length scale. Now 
looking back at the pair correlation function (Figj^), the correlation peaks beyond 10 nm 
scale should arise from the inter-cluster interaction. No matter how strong the lamellar peak 
is, the mass fractal accompanied by clustering is an universal property in spider dragline 
silk, manifested by the power law of structure factor and the reconstructed electron density 
map (see Fig.S3). 

The mass fractal property is signihcant to the mechanical strength of spider silks. The 
clusters can be viewed as the basic functional mechanical units, and more importantly 
such functional units exist throughout the intermediate length scale. Whenever the silk 
experiencing an external force, the largest cluster of the size 100 nm will deform to respond to 
the external kick, then the deformation is propagated down to the smaller clusters which will 
deform in response. The force-induced deformation will be propagated down in this manner 
all the way to the /S-sheets level, causing an exponentially increasing number of deformed 
functional units as the length scale reduces. As illustrated by Figj^, the hierarchical self¬ 
similar model is remarkably robust. This scheme is consistent with the model proposed by 
Zhou & Zhang [21] which effectively explained the exponential force-extension property of 
spider silks To answer the questions imposed by the large crystalline regions observed 
in the SEM and TEM images, we suggest that these 20 - 50 nm granules are composed of 
densely packed, uniformly oriented and strongly interconnected /3-sheet crystals of the sizes 
2 to 4 nm. 
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III. EXPERIMENTAL METHODS 


Major ampullate dragline silks were collected by forced silking from living spiders anes¬ 
thetized with carbon monoxide. The silks were reeled at 2.0 ± 0.1 cm-s'^ and were directly 
mounted across hollow cardboard holders driven by an electric motor. At these silking 
speeds, /3-sheet crystals are supposed to be extremely well oriented |23]- The samples nor¬ 
mally have 50 to 100 strands of hbres. 

The X-ray experiment was carried out at Argonne National Laboratory Advanced Photon 
Source BioCars 14-ID-B beam line optimized with a SAXS setup. The incident beam energy 
was 9 KeV, corresponding to an X-ray wavelength of 1.38 A. The scattering data were 
collected on a Marl65, a 3072 by 3072 pixel resolution CCD detector from Marresearch with 
a pixel size of 79 pm. The detector to sample distance is hxed at 180 mm throughout the 
experiment. The sample was mounted to an xyz-translation goniometer head with helium 
gas chamber placed between the CCD detector and the sample to help reduce the background 
signal. The exposure time varied from 2 s to 6 s depending on the size of the silk bundle. 
For each sample ten exposures were taken. Backgrounds were measured after translating 
the silk bundle out of the X-ray beam. 

The multiple X-ray exposures were averaged and background subtracted by using software 
Fit2D [21]. The integrated 1-D WAXS and SAXS prohle were obtained by using azimuthal 
integration functionality in Fit2D. The automatic Gaussian peak htting was performed using 
the Isqcurvefit optimization routine from MATLAB. 


IV. SIMULATION METHOD 

The transformation from electron density map to scattering intensity adopts the fast 
Fourier transform (FFT) simulation method proposed by Klauss Schmidt-Rohr [T2]. The 
electron density map p(r) was represented by an V x V matrix, where N is usually chosen 
to be a power of 2, and then converted to a reciprocal space scattering intensity /(q) by 
the transformation algorithm. For each sample, we initialized the silk structure such that 
it conforms to the crystal sizes calculated from the WAXS data and maintained the closest 
approach distances, which is the d-spacing of SAXS lamellar peak (Tab.l). The scattering 
centers, which are the /3-sheet crystals in this model, were represented by higher contrast 
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rectangular shape sub-matrix and are uniformly oriented to be parallel to the silk hber axis. 
The orientation of the crystals at >2.0 cm-s'^ can be approximated by perfect alignment 
with very small error [2T1132]- A hard-shell exclusion geometry was used to maintain the 
closest approach distance. 

We imposed a model with a lamellar modulation by generating 60 nm wide crystal-rich 
stripes separated by equal width empty spaces. The lamellar structure is essential to physi¬ 
cally represent the silk hbrils hne structure within the silk hber [T2ll23]-The initial crystalline 
map was then fed to the stimulated annealing reconstruction routine. The reconstruction 
was proceeded by generating a electron density map such that the calculated structure factor 
S{q) matches the experimental S{q) with acceptable error tolerance. This was achieved by 
minimizing the pseudo-energy 


Er.e^ = '£\S{<l)-S(q)\'‘ (4) 

Q 

which measures the distance from simulated structure factor to the experimental value [Mh 
155] . The optimization was realized by the stimulated annealing algorithm where each random 
walk is accepted or rejected by a probability of [5B] 

P(s' s) = (5) 

ksT 

where s and s' are the states before and after one random walk, AEgi^g is the change 
of pseudo-energy after accepted state transition, ks is the Boltzmann constant, T is the 
imaginary stimulated annealing temperature. 


V. ACKNOWLEGMENTS 

This work was supported by the Department of Defense, AFOSR (FA9550-14-1-0014) and 
the US National Science Foundation (DMR-1264801). The author would like to thank Klaus 
Schmidt-Rohr for providing the MATLAB simulation code, Robert Henning at APS BioCars 
for setting up 14-ID-B beamline and Yang Jiao for instructive discussion. This research used 
resources of the Advanced Photon Source, a U.S. Department of Energy (DOE) Office of 
Science User Facility operated for the DOE Office of Science by Argonne National Laboratory 
under Contract No. DE-AC02-06CH11357. Use of the BioCARS Sector 14 was supported 


13 



by the National Institutes of Health, National Center for Research Resources, under grant 
number RR007707. 


[1] L. Romer and T. Scheibel, Prion 2, 154 (2008). 

[2] R. V. Lewis, Chemical Reviews 106, 3762 (2006), 

[3] S. Keten, Z. Xu, B. Ihle, and M. J. Buehler, Nat Mater 9, 359 (2010). 

[4] J. E. Jenkins, S. Sampath, E. Butler, J. Kim, R. W. Henning, G. P. Holland, and J. L. Yarger, 
Biomacromolecules 14, 3472 (2013). 

[5] T. Asakura, Y. Suzuki, Y. Nakazawa, G. P. Holland, and J. L. Yarger, Soft Matter 9, 11440 
(2013), 

[6] C. Y. Hayashi, N. H. Shipley, and R. V. Lewis, International Journal of Biological Macro¬ 
molecules 24, 271 (1999). 

[7] M. Xu and R. V. Lewis, Proceedings of the National Academy of Sciences 87, 7120 (1990), 

[8] S. Keten and M. J. Buehler, Phys. Rev. Lett. 100, 198301 (2008), 

[9] S. Keten and M. J. Buehler, Journal of the Royal Society Interface 7, 1709 (2010). 

[10] G. P. Holland, Q. Mou, and J. L. Yarger, Chem. Commun. 49, 6680 (2013), 

[11] C. Riekel and E. Vollrath, International Journal of Biological Macromolecules 29, 203 (2001). 

[12] S. Sampath, T. Isdebski, J. E. Jenkins, J. V. Ayon, R. W. Henning, J. P. R. O. Orgel, 
O. Antipoa, and J. L. Yarger, Soft Matter 8, 6713 (2012), 

[13] M. Kitagawa and T. Kitayama, Journal of Materials Science 32, 2005 (1997), 

[14] J. E. Trancik, J. T. Czernuszka, C. Merriman, and C. Viney, Journal of Microscopy 203, 235 

( 2001 ), 

[15] A. L. Patterson, Phys. Rev. 56, 978 (1939), 

[16] D. J. Yarusso and S. L. Cooper, Macromolecules 16, 1871 (1983), 

[17] K. Schmidt-Rohr and Q. Chen, Nat Mater 7, 75 (2008). 

[18] K. Schmidt-Rohr, Journal of Applied Crystallography 40, 16 (2007), 

[19] J. S. Pedersen, Journal of Applied Crystallography 27, 595 (1994), 

[20] M. S. Creager, J. E. Jenkins, L. A. Thagard-Yeaman, A. E. Brooks, J. A. Jones, R. V. Lewis, 
G. P. Holland, and J. L. Yarger, Biomacromolecules 11, 2039 (2010), 


14 














[21] M. Lei, A. M. R. de Graff, M. F. Thorpe, S. A. Wells, and A. Sartbaeva, Phys. Rev. B 80 , 
024118 (2009), 

[22] T. Proffen and S. J. L. Billinge, Journal of Applied Crystallography 32 , 572 (1999). 

[23] N. Du, X. Y. Liu, J. Narayanan, L. Li, M. L. M. Lim, and D. Li, Biophysical Journal 91, 
4528 (2006). 

[24] K. J. Koski, P. Akhenblit, K. McKiernan, and J. L. Yarger, Nat Mater 12 , 262 (2013). 

[25] J. Teixeira, Journal of Applied Crystallography 21 , 781 (1988), 

[26] J. E. Martin and A. J. Hurd, 

[27] D. W. Schaefer, Science 243 , 

[28] H. Stanley, Journal of Statistical Physics 36, 843 (1984). 

[29] H. Zhou and Y. Zhang, Phys. Rev. Lett. 94, 028104 (2005). 

[30] N. Becker, E. Oroudjev, S. Mutz, J. P. Cleveland, P. K. Hansma, C. Y. Hayashi, D. E. 
Makarov, and H. G. Hansma, Nat Mater 2 , 278 (2003). 

[31] A. P. Hammersley, S. O. Svensson, M. Hanfland, A. N. Pitch, and D. Hausermann, High 
Pressure Research 14 , 235 (1996), 

[32] J. D. van Beek, S. Hess, E. Vollrath, and B. H. Meier, Proceedings of the National Academy 
of Sciences 99, 10266 (2002), 

[33] R. L. McCreevy and L. Pusztai, Molecular Simulation 1 , 359 (1988), 

[34] R. Kaplow, T. A. Rowe, and B. L. Averbach, Phys. Rev. 168 , 1068 (1968), 

[35] Y. Jiao, E. H. Stillinger, and S. Torquato, Proceedings of the National Academy of Sciences 
106 , 17634 (2009), 

[36] S. Kirkpatrick, C. D. Gelatt, and M. P. Vecchi, Science 220 , 671 (1983), 


Journal of Applied Crystallography 20 , 61 (1987) 
pp. 1023 (1989). 


15 














